5 



X 



UK/93-02 

Origin of Difference Between d and u Partons in the 

Nucleon 



m 

0\ ■ Keh-Fei Liu and Shao-Jing Dong 



Department of Physics and Astronomy 
^ ! University of Kentucky 

Lexington, KY 40506 



(N 

J> I Abstract 

0^ 

0^ ■ Using the Euclidean path- integral formulation for the hadronic tensor, we 

. show that the violation of the Gottfried sum rule does not come from the 

I disconnected quark-loop insertion. Rather, it comes from the connected (quark 

' line) insertion involving quarks propagating in the backward time direction. 

We demonstrate this by studying sum rules in terms of the scalar and axial- 
vector matrix elements in lattice gauge calculations. The effects of eliminating 
backward time propagation are presented. 



g^' PACS numbers: 13.60.Hb, ll.15.Ha, 12.38.Gc 

A recent measurement of the Gottfried sum rule(GSR), defined in terms of the 
difference between the proton and neutron structure functions F2{x) in the integral 
Sc = /o dx[Fi{x) - F^{x)]/x, by the New Muon Collaboration (NMC) has shown 
5^ ! a disagreement with the expectation of the naive parton model. Assuming charge or 

isospin symmetry, the sum rule Sq can be expressed in terms of the parton distribu- 
tions in the parton model as 



5g = ^ + ^ ['dx[uP{x)-dPix)]. 
6 6 Jo 



The naive parton model which assumes the isospin symmetry in the "sea", i.e. 
u'^{x) = d^{x), leads to the prediction that So = 1/3 0. However, the NMC data, 
extrapolated to a; = and 1, leads to a value So = 0.24 ± 0.016 which implies that 
vP{x) and d^{x) are not the same in the proton with the number of vP less than that 
oidP. 

This has generated a good deal of theoretical interest. The apparent isospin 
asymmetry in the sea was envisioned by Field and Feynman P] as due to the Pauli 
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exclusion principle and has been modeled Q with the Sullivan process which 
considers the meson cloud in the nucleon and the chiral-quark model 0. 

In order to gain insight into the origin of this large d/u difference in the "sea", we 
will examine the deep inelastic scattering in the Euclidean path integral formalism. 
The advantage of this formalism is that one can follow the quark line of propagation in 
the Euclidean time and separate out different contributions in terms of the connected 
and disconnected quark line insertions to facilitate the discussion. 

The deep inelastic scattering of muon on nucleon involves the hadronic tensor of 
the current- current correlation function in the nucleon, i.e. 

W^u{q^,v) = j^{N\ J ^e''^-''J^{x)J^{0)\N),pi„ave.- This forward Compton ampli- 
tude can be obtained by considering the ratio of the four-point function 
{ONit)J^ix,ti)J^{0,t2)ONiO)) and the two point function {ONit- {h - t2))C>Jv(0)), 
where ON{t) is the interpolation field for the nucleon at Euclidean time t with zero 
momentum. For example, Oi\i{t) can be taken to be the 3 quark fields with the nucleon 
quantum numbers. On = J d'^xe"''"^'^^^^"-{x){{'^^^^''{x))'^C'y5'^^'^^'^{x)) for the proton. 

As both t — ti >> I/AMtv and t2 >> I/AMtv, where AM^ is the mass gap 
between the nucleon and the next excitation (i.e. the threshold of a nucleon and a 
pion in the p-wave), the intermediate state contribution in the four-point and two- 
point functions will be dominated by the nucleon with the Euclidean propagator 
g_Mjv(t-(ti-t2))_ Hence, 



< O(t)/|fe-^f%(f,ti)J,(0,t2)O(0) > 



< O{t-T)O{0) > 



t-ti» 1/AMn 
t2» 1/AMn 



2Mn 



p^-MM{t-r) 

1 < iV| / ^e-^^"%(x,ti)J,(0,t2)|iV >, 
J Ztt 
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where r = ti — f is the transition matrix element (O|07v|-^), and V is the 3- volume. 
The hadronic tensor can be obtained formally by the inverse Laplace transform [0, 
W^uiq"^, ^) = j Ic-i^ dre'^^W^uiq'^, r) or through the integration 
W^Aq'', ^) = |lim,^oi?e/o'£r2e(^+^^)^l^^,(g2, r)rfr with c > 0. 

In the Euclidean path-integral formulation, the four-point function can be classi- 
fied into different groups depending on different topology of the quark paths between 
the source and the sink of the proton. They represent different ways the fields in the 
currents and Jy contract with those in the nucleon interpolation operator Oat at 
different times. This is so because the quark action and the electromagnetic currents 
are both bilinear in quark fields, i.e. in the form of \E'M\E', so that the quark numbers 
are conserved and as a result the quark line does not branch the way a gluon line 
does. As illustrated in Fig. 1, we see Fig. 1(a) and 1(b) represent connected inser- 
tions (C.I.) where the quark fields from the currents contract with those from On 
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and the quark lines from t = to t = t are connected with the currents. Fig. 1(c), 
on the other hand, represents a disconnected insertion (D.I.) where the quark fields 
from Jn and Ji, self-contract and are hence disconnected from the quark paths which 
originate from On{0) and terminate at Oiv(t). Here, "disconnected " refers only to 
the quark lines. Of course, quarks sail in the background of the gauge field and all 
quark paths are ultimately connected through the gluon lines. The infinitely many 
possible gluon lines and additional quark loops are implicitly there in Fig. 1 but are 
not explicitly drawn. Fig. 1 represent the contributions of the class of "handbag" 
diagrams where the two currents are hooked on the same quark line. These are lead- 
ing twist contributions in deep inelastic scattering. The other contractions involving 
the two currents hooking on different quark lines are represented in Fig. 2. Given a 
renormalization scale, these are higher twist contributions in the Bjorken limit. We 
shall neglect these "cat's ears" diagrams in the following discussion. 

In the deep inelastic limit where < 0(1/(5^) (we are using the Minkowski nota- 
tion here), the leading light-cone singularity of the current product (or commutator) 
gives rise to free quark propagator between the currents. In the time-ordered dia- 
grams in Fig.l, Fig. l(a)/l(b) involves only quark/antiquark propagator between 
the currents. Whereas, Fig. 1(c) has both quark and antiquark propagators. Hence, 
there are two distinct classes of diagrams where the antiquarks contribute. One comes 
from the D.I. ; the other comes from the C.I.. It is frequently assumed that connected 
insertions involve only "valence" quarks which are responsible for the baryon number. 
But apparently, this is not true. To define the quark distribution functions more pre- 
cisely, we shall call the antiquark distribution from the D.I., which are connected to 
the other quark lines through gluons, the sea antiquarks and the antiquarks from the 
C.I. the cloud antiquarks [Q. Thus, in the parton model, the antiquark distribution 
function can be written as 

q'ix)=ql{^)+t{^). (3) 

to indicate their respective origins in Fig. 1(b) and Fig. 1(c) for each fiavor i. 
Similarly, the quark distribution will be written as 

q\x) = qUx) + qlix) + qlix) (4) 

where ql{x) comes from Fig. 1(c) and both qy{x) and ql{x) come from Fig. 1(a). Since 
ql{x) = q^si.x)i define q],{x) = q^^i^x) so that qy(x) will be responsible for the baryon 
number, i.e. / uv{x)dx = J[u{x) — u{x)]dx = 2 and / dv{x)dx = J[d{x) — d{x)] = 1 
for the proton. 

We shall first examine Fig. 1(c). After the integration of the Grassman fields \E' 
and \E', the path-integral for Fig. 1(c) can be written as the correlated part of 

J rf[A]e-^«Tr[M-i(t2,ti)7.M-i(ti,t2)7M]r^[^"'(^>0)...M-i(t,0)...M-i(t,0)...]. 

(5) 

where A is the gluon field. So the gluon action, and M is the quark matrix in the 
bilinear quark action \&M\E'. M~^{ti,t2) denotes the quark propagator from t2 to ti. 
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Note in eq.(^, the trace is over the color-spin as well as the flavor indices. Since the 
quark loop involving the currents is separately traced from those quark propagators 
M~^(t,0) whose trace reflects the quantum numbers of the proton, eq.(^ does not 
distinguish a loop with the u quark from that with the d quark at the flavor-symmetric 
limit, i.e. m„ = rrid. These are referred to as sea quarks and sea antiquarks in the 
naive parton model, since they are connected to those quark propagators which are 
sensitive to the hadron state through the gluon lines. These sea quarks can not give 
rise to the violation of the GSR, since Ug = dg- The isospin breaking will give a small 
effect in the order of {rriu — m(i)/Mc 0, where is the constitute quark mass which 
reflects the confinement scale. Hence, the isospin symmetry breaking effect will be 
at the 1% level. It does not explain the violation of the GSR which is at ~ 30% 
level On the other hand, the quark propagators connecting the currents in Fig. 
1(b) will show up in the same trace along with other quark propagators connecting 
the interpolation fields. Therefore, the cloud antiquarks are subjected to the Pauli 
exclusion as does the valence quarks and cloud quarks in Fig. 1(a) [jlO[. Consider 
the Fock space where a u quark line does the twisting in Fig. 1(b), the simplest 
Fock space would then be uuuud. With 3 u quarks, this Fock space configuration 
might be more Pauli suppressed than the corresponding Fock space of uuddd with 2 
u quarks and 2 d quarks. We believe this is the reason for the large d/u difference 
in the nucleon as revealed by the NMC data. Consequently, neglecting the isospin 
symmetry breaking, the sum rule 5*^ can be written as 

1 2 /"i - 
Sg = - + - / dx^c{x) - d^{x)\ (6) 
3 3 Jo 

How do we substantiate this claim? Instead of evaluating the hadronic tensor 
directly which involves a four-point function, we shall study matrix elements with one 
current which can be obtained from three point functions. In the spirit of the operator 
product expansion and the parton model, matrix elements of the twist-2 operators 
in the form (A^lvE^PvE'lA^) are the sum rules of the parton distribution functions. This 
can be viewed as — >■ in the Bjorken limit, the two currents at ti and t2 merge 
into one so that the connected insertion of one local operator will have both types 
of paths represented in Fig. 1(a) and Fig. 1(b). In fact, there have been indirect 
evidences of the presence of the cloud antiquarks in the previous study of three point 
functions in the quenched lattice QCD calculations, such as the p meson dominance 
in the pion electric form factor ||ll|] and the negative neutron charge radius |jl2[ . It has 
also been considered in association with large N^. and chiral perturbation theory 



To explicitly reveal the existence of the cloud antiquarks, we shall consider the scalar 
and axial-vector matrix elements in lattice calculation. The scalar current expanded 
in the plane-wave basis 

d'^x^^ix) = J d^k— J2[bHk, s)b{k, s) + d\k, s)d{k, s)] (7) 



4 



is a measure of the quark and antiquark number up to the factor m/E. Since the 
first moment of the structure function F2 is not expressible in terms of the forward 
matrix element of a leading twist-2 operator, the scalar matrix element has been 
taken as a measure of the quark and antiquark number with m/E approximated 
by a constant [0, |T^. For the lattice calculation, we shall consider quark masses 



ranging from the charm to strange. In this case, we expect the dispersion of E 
will be considerably narrow so that m/E will be close to unity. To further decrease 
the dependence on m/E and other lattice corrections like the finite volume effect, 
scaling, and finite lattice renormalization, we shall consider ratios of matrix elements. 
Furthermore, since we have shown that the sea quarks from the D.I. can not give 
any significant contribution to the GSR, we shall concentrate on the C.I.. The first 
ratio we calculate is the isoscalar to isovector forward scalar matrix element or scalar 
charge of the proton with C.I.. In the parton model, it should be written according 
to eqs. (13) and (H) as 



^ {p\uu\p) - {p\dd\p) 



{p\uu\p) + {p\dd\p) 



1 + 2 / dx[uc{x) — dc(x)] 
C.I. ~ 3 + 2 J dx[ucix) + d^{x)] 



In view of the fact that 5*^ shows / dx[uc — dc] < experimentally, we expect this 
ratio to be < 1/3. Our lattice results based on quenched 16^ x 24 lattice with j3 = 6 
for the Wilson k ranging between 0.154 to 0.105 which correspond to strange and 
twice the charm masses are plotted in Fig. 3. For heavy quarks (i.e. k > 0.140 
or rUga > 0.31 in Fig. 3), the ratio is 1/3. This is to be expected because the 
cloud antiquarks which involves Z-graphs are suppressed for non-relativistic quarks 
by 0{p/mq). For quarks lighter than /t = 0.140, we find that the ratio is in fact 
less than 1/3. We take this to be the evidence of the cloud antiquarks in eq. (||). 
To verify the fact that this is indeed caused by the backward time propagators, we 
perform the following simulation. In the Wilson lattice action, the backward time 
hopping is prescribed by the term —k{1 — '~^4)Ui{x)5x,y-a4,- We shall amputate this 
term from the quark matrix in our calculation of the quark propagator. As a result, 
the quarks are limited to propagating forward in time and there will be no Z-graph 
and hence no cloud quarks and antiquarks. The Fock space is limited to 3 valence 
quarks. This is what the naive quark model is supposed to describe by design. In this 
case, the scalar current in eq. and the ratio in eq. (|^) involve only the valence 
quarks. To the extent the factor m/E can be approximated by a constant factor, the 
ratio Rs in eq. (|^) should be 1/3. The lattice results of truncating the backward time 
hopping for the light quarks with k = 0.148, 0.152 and 0.154 are shown as the dots in 
Fig. 3 with errors less than the size of the dots. We see that they are indeed equal to 
1/3. This shows that the deviation of Rs from 1/3 is caused by the cloud quarks and 
antiquarks. In retrospect, this can also be used to justify approximating m/E by a 
constant factor in eq. (^. We also find that the isovector scalar charge of the proton 
(the numerator in eq. (H)) for the forward propagating case is greater than the case 
with both forward and backward time propagation in our lattice results. For instance. 
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the isovector scalar charges for the forward propagating case are 1.07(1) and 1.01(1) 
for K = 0.154 and 0.152. Yet, they are 0.73(15) and 0.85(5) respectively for the case 
with both forward and backward propagations. Assuming the m/E factor to be the 
same for these two cases and other things being equal, it implies that / dx[uc — dc] < 
which is consistent with the NMC result. 

We have also examined the ratio of the isoscalar to isovector axial charge (C.I. 
only) of the proton. In the parton model, the ratio can be written as 



Ra 



{p\u'y3l5u\p) + {p\d^3l^d\p) g\ Jdx[Au{x) + Ad{x)] 



{p\u-f3l5u\p) - {p\d-f375d\p) 



C.I. g\ 



C.I. S dx[Au{x) - M{x)] 



C.I. 

(9) 

where Au(A(i) is the quark spin content of the u(d) quark and antiquark in the C.I. 
At the non-relativistic limit, is 5/3 and g\ for the C.I. is 1 (the spin of the proton 
is entirely carried by the quarks in this case) |T6[. Thus, the ratio should be 3/5 and 
we did find this ratio for the heavy quarks in Fig. 4. For lighter quarks, the ratio 
dips under 3/5. We interpret this to be due to the cloud quark and antiquark. Again 
when we dropped the backward time propagation for the quarks, we find that the 
ratio shown as the dots in Fig. 4 becomes 3/5 for lighter quarks as predicted by the 
quark model. 

There are phenomenological consequences for the cloud quarks and antiquarks. 
The structure functions extracted from the DIS need to reflect the definitions of the 
quark and antiquark distributions in eqs. @) and (^). Since strange and charm quarks 
come only from the sea in Fig. 1(c), it is natural to expect that u,d > s,c since the u 
and d have both the sea and the cloud parts. The neutron-proton mass difference can 



be understood in terms of the isovector scalar charge in eq. The violation 
of the GSR has been modeled in terms of the Sullivan process 0| and the chiral 
quark model Although these models give the right picture in terms of the cloud 
antiquarks, there are inevitable drawbacks in the effective theories. For example, the 
Sullivan process where the photon couples to the antiquark in the meson as depicted 
in Fig. 5(a) can be drawn in terms of the quark lines in Fig. 5(b). However, Fig. 
5(b) is only half of the story as far as the forward Compton amplitude is concerned. 
Upon attempting to complete the other half, one has a choice of taking the mirror 
image of Fig. 5(b) which will lead to the D.I. in Fig. 1(c) which does not contibute 
to the GSR. Alternatively, one could sew the quark lines with one of them twisted 
which will then lead to Fig. 1(b) and this does contribute to the GSR. The Sullivan 
process and the chiral quark model do not distinguish these two different topological 
possibilities. 

In conclusion, we have shown in the Euclidean path-integral formalism that the 
experimentally observed d/u difference in the proton comes from the C.I. which in- 
volves cloud quarks and antiquarks. We have studied it in terms of the ratios of 
the isovector to isoscalar scalar and axial charges of the proton in the lattice calcu- 
lations for the C.I.. We found that these ratios have the expected non-relativistic 
and relativistic limits as far as the cloud antiquarks are concerned. We demonstrate 
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this by truncating the quark backward time propagation which leads to the quark 
model predictions for these ratios without cloud antiquarks. Other phcnomcnological 
implications related to the cloud quarks and antiquarks and the quark model will be 
explored in the future. 
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Figure Captions 

Fig. 1 Time-ordered "handbag" skeleton diagrams of quark lines with different topolo- 
gies, (a) / (b) is the C.I. involving a quark/antiquark propagator between the currents, 
(c) is a D.I. involving sea quarks and antiquarks. 
Fiq.2 Cat's cars diagrams. 

Fig. 3 The ratio of the isovector to isoscalar scalar charge of the proton for the C.I. 
(shown as O) is plotted as a function of the quark mass ruq in the lattice unit a. The 
errors are obtained from the jackknife method. The errors of the dots are smaller 

than the size of the dots. 

Fig. 4 The ratio of the isoscalar to isovector of the proton for the C.I. as a function 
of the quark mass. 

Fig. 5 (a) Sullivan process in terms of the meson and baryon lines, (b) The same 
process drawn in terms of the quark lines. 
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